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Abstract: Considering the equality of the proton and neutron Fermi levels as 
^vq ' an indication of a phenomenological interaction between pairs of protons and neu- 

£jry trons, low-energy isovector monopole resonances are proved to appear in the super- 

fluid nuclei. The required phenomenological interaction is presented as an isospin 
symmetry-breaking mean field for the four particle interaction. 

in 



1. Introduction 

> 

' The isovector monopole resonances (IVMR) were predicted in the frame- 

work of the nuclear hydrodynamic model ^ [2] as out-of-phase breathing 
, oscillations of the proton and neutron densities. These oscillations repre- 

sent the isovector counterpart to the isoscalar compression modes, and are 
generated by the restoring force which appears from the symmetry energy. 
Experimentally, IVMR were first observed in the charge exchange reactions 
90 Zr, 120 Sn (tt~ , ir°) at T v - = 165 MeV [3j, and were also confirmed recently 
in the (n,p) reactions on 90 Zr [I]. Other experiments J5] have shown their 
occurrence for a wide range of nuclei, from 40 Ca to 208 Pb, both normal and 
superfluid. 

The quantum description of the superfluid systems shows a close con- 
nection between the velocity potential x, (v = — Vx); describing the flow 
and the gauge angle ip from the BCS transformation to quasiparticles [2]: 
V(x — ip/2m) = 0, where m is the nucleon mass. For the monopole vi- 
brations, this connection implies a radial dependence of the proton and 
neutron gauge angles (p p , (p n . In contrast to the space dependence, the time 
dependence of tp T , r = p,n occurs even if no macroscopic flow is present 
(Vx = 0). The time-derivative (p T /2 is just the Fermi energy A T and a 
superfluid system performs in its ground state a free gauge rotation with the 
angular velocity <p T = 2A r . As is known, the proton and neutron systems 
are not closed, changing particles through f3 decay up to the equalisation 
of the Fermi energies A p , A n . This can be taken as an indication of a phe- 

1 first printed as preprint FT-353-(1989)/April, by the Institute of Physics and Nuclear 
Engineering, Bucharest, Romania. 
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nomenological gauge-restoring interaction, tending to fix not only the space- 
depending term, but the whole relative angle (p p — ip n at a constant value. 
Consequently, low-lying IVMR, appearing only in superfluid nuclei are ex- 
pected, as the isovector counterpart to the isoscalar uniform gauge rotation. 
The phenomenological character of the gauge-restoring interaction gener- 
ating these low-lying modes is imposed by its non-commutation with the 
particle-number operators N p , N n . So, it can be considered only as a mean- 
field approximation for a microscopic number-conserving Hamiltonian. 

In this paper section 2 presents the phenomenological gauge interaction 
as a mean field deriving from a microscopic four-particle interaction. Restor- 
ing the isospin symmetry by the cranking method, a kinetic energy term 
contributing to the symmetry energy is obtained. These connections be- 
tween the gauge-restoring interaction, the four-particle interaction, and the 
symmetry energy are used at the end of section 2 to estimate the former's 
strength 5$. Following the semiclassical treatment of the proton- neutron 
quadrupole-quadrupole interaction at deformed nuclei [Hj in section 3 the 
effects of the gauge interaction in superfluid nuclei are investigated. This 
approach gives both the underlying classical picture and the quantised spec- 
trum of the Fermi levels oscillations, allowing for realistic estimates of their 
energy. 

2. The gauge-restoring interaction 

The early attempts to construct an isoscalar four-particle interaction 
were based on the algebraic properties of the pairing operators, of hav- 
ing closed commutation relations with the isospin operators. Denoting by 
(a,m), a = (n,l,j), the shell-model states (n,l,j,m), the proton-proton, 
neutron- neutron and the proton-neutron pair creation operators P+, N + , 
R + , are: 






(3) 



If P_, P_ are their hermitian conjugates, then 



[P+,P_] = 2P 



[JV+,J\L] =2N 



[R+, R-] — 2Rq 



(4) 
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Figure 1. Root diagram of o(5) algebra, 
with 

^0 = o 'y^^ C r>arn c pam ~ r) (5) 
a,m 

-^0 Q ^ ] ( c vam c 'parn + Cnam c nam ~ 1) (7) 

each couple (P_, Po, P+), (JV_, JVo,iV+), (P-, Po, P+) generating an su(2) 
algebra. Moreover, if (T_,To,T + ) are the isospin operators: 

P+ = c Lm c P«m ) T— = c pam c nam (8) 
a,m a, m 

-^0 = 2 '^ / ( C nam c nam ~ c pam c pam) (9) 

then (P±, iV±, P±, T±, Po, iVp) generate the second-rank semisimple Lie alge- 
bra o(5) [7j. Choosing Po and To as basis elements in the Cartan subalgebra, 
the root diagram can be pictured as in figure 1 |S]. 

The three commuting operators V_ x = P+, V\ = R+/y/2, v\ = N + 
are the components of an isospin vector, and were used in |3 to con- 
struct the quadrupling four-particle interaction Hq = —GqQqQq/A, with 
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Qo = ^E^it"!)' 1 ^-/!' This interaction does not lead to the expected 
gauge force for the superfiuid systems, and consequently in the present work 
a different choice will be made. Coupling the isovector with its Hermi- 
tian conjugate to an isospin quadrupole Qi^-, an isoscalar four-particle 
interaction, separable in the two particles-two holes channel will be defined 
by 

2 2 

H 4 = -G 4 E (-1) M Q2 M Q2-^ = -G? 4 E QlnQtn ■ ( 10 ) 

[m=-2 n=-2 

With this term, the total microscopic Hamiltonian, including the usual pair- 
ing interaction, becomes 

H = H + H 4 (11) 

-^0 = ^ A ^pgCyamCpam ~l~ ^naCnam^nam) GpP+P— G n N+N— . (12) 

a,m 

The ground state \g) of H includes the 4-particle correlations determined 
by H 4 and the transition to the gauge-restoring mean field is expected when 
the sum rule: 

2 2 

E (9\QlQ2,\g)= E ElH^b)l 2 (13) 

H=-2 n=-2 \ n ) 

is exhausted by |^): 

2 2 

E (9\Ql^\g) « E I?/ > ( 14 ) 

//=-2 //=-2 

^ = (g\Q2M (15) 

and 5-J-2 7^ 0. This highly correlated ground state can be approximated by 
the eigenfunction \gj) of the linearised Hamiltonian H^: 

2 

H L = H - G 4 E Ml + QW ~ " T o (16) 

H=-2 

H L \g w ) = E u \g u ) (17) 
where are self-consistently determined by: 

% = (gw\Q2»\gJ) ■ (18) 

The "cranking" term —loTq was introduced for an approximate projection of 
the proton and neutron numbers Z, N, because even if Hi commutes with 
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A = N p + N n , it does not commute with To. So, the parameter to is fixed 
by the constraint: 

N - Z 

{gu\T \g u ) = —j— ■ (19) 

To obtain the function \gu) further approximations will be made, neglecting 
the proton-neutron pairing interaction and retaining in \gj) only the corre- 
lations between pairs of protons and pairs of neutrons. Consequently, Hl 
becomes: 

H L = ^a4 am c ram - G P P+P_ - G n N+N- + H g - ljT , (20) 

r,a,m 

where 

G T = G T + (21) 

and H g = -2G 4 (q%N + P_ + q 2 N^P + ), q 2 = {g w \N + P„\ gu> ) is the gauge 
restoring interaction. This Hamiltonian will be restricted to one degenerate 
j = 3/2 level. The interest for this case is justified by the relevance of 
the 0(5) symmetry for the classification of the experimental data. So, the 
low- lying + states of the nuclei filling j = 3/2 shells, as IP3/2, ld3/2 and 
2P3/2 can b e assigned to the weight vectors of the 14-dimensional irreducible 
representation of the o(5) algebra. In these multiplets both ground and 
excited states are included, and figure 2 presents their relative energies with 
Coulomb corrections 2 . Considering the same pairing constant G and single- 
particle energy e for both protons and neutrons, the restricted Hamiltonian 
is 

H L = 20e + (2e + co)P + (2e - w)N (22) 

-G(P + P^+N + N_)-2G 4 (q* 2 N + P-+q2N-P + ) , ^ = 3 + \ ■ 

This Hamiltonian commutes with the operator A and with the squares of the 
proton and neutron quasispin operators P 2 , iV 2 , so that its eigenfunctions 



2 The experimental binding and excitation energies have been taken from A. H. Wapstra 
and K. Bos, The 1977 atomic mass evaluation, Atomic Data and Nuclear Data Tables 19 
177 (1977), and P. M. Endt and C. van der Leun, Energy levels of A=21~44 nuclei, Nucl. 
Phys. A 214 1 (1973). When I obtained figure 2 (after the less symmetric diagram for the 
light nuclei filling the IP3/2 shell), I wondered if there was no previous attempt to make 
such a test of the 0(5) symmetry. After a thorough search in literature for several weeks 
I found a practically unknown paper, by P. Camiz and U. Catani |11|. in which the 0(5) 
symmetry was used to classify nuclear states for a wide range of mass numbers. However, 
I could not find any indication that an effective proton-neutron interaction like H g of 1201 1 
has ever been used before. 
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Figure 2. The o(5) multiplet of J ff = + ground (solid) and excited (dash) states for 
the nuclei filling the ld 3 / 2 valence level. Whenever available, the excitation energy 
with respect to the ground state ((2+; 1) for 36 C1, 36 K and (3+; 0) for 38 K) is given 
in round brackets (MeV). 



are labelled by the eigenvalues E, A v = number of the valence particles, p 
and n: 

H L \pnEA v )u = Ej\pnEA v )u (23) 
A\pnEA v )u = A v \pnEA v )u (24) 
P 2 \pnEA v ) w = p(p+l)\pnEA v )u (25) 
N 2 \pnEA v ) w = n(n + l)\pnEA v )u . (26) 

In addition, {pnEAv}^ depends on oj whose value is fixed by the constraint 
(19): 

TV — Z 

u {jmEA v \T Q \pnEA v ) u = v ^ v . (27) 

So, the function IpnEAv)^ is associated with a given nucleus from the iso- 
baric family determined by A v . The diagram in figure 2 contains five such 
families, having in the increasing order of the energies: A v = 8, 6, 4, 2, 0. 
From these, the correlation between the proton and neutron pairs is effective 
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only for p = n = 0/2 = 1, A v = 20 = 4, and this case will be discussed in 
the following. Denoting by \p,po) and \n,no) the seniority eigenfunctions of 
P 2 , Pq and N 2 , Nq, respectively, the function IpnEAv)^ can be expressed 
as: 

\pnEA v ) LU = a po \p,p Q )\n,n ) . (28) 

Po+no=A v /2— fi 

When p = n = 1, A v = 4, = 2, the coefficients a_i, ao, oi, are easily 
found from the equation iTkllliM)^ = £ , a ,|ll£ , 4} £J , and (|28f) becomes 



| 1 1 t I ) ,. 1 : ( -^L_ |! I .1 (29) 



-|1,0>|1,0> + 1>|1,-1» . 

Here x = 2G-\-E LlJ and 5^ = 2G*4q , 2, with (72 chosen to be real, self-consistently 
determined by 

8x5 

q 2 = w (lLE4|iV +J P_|ll£4) w = j- A a . (30) 

(^-^(l + S^^^) 

The w-dependent energy £7 W is analitically obtained from the third-order 
equation 

x 3 + 2Gx 2 - 4{25l + oj 2 )x - 8Guj 2 = . (31) 

Its solutions are 

E *>=-^--\ cos yG 2 + 65 2 + 3u; 2 , (32) 

E ku =~ + |(«>s I + (-l) fc V3sin ^& + 65 2 + 3u; 2 , fc = l,2 (33) 
with ^ G [0, 7r] , 

tan<£ = - , (34) 

7 



/fi4 r* 2 

« = \/^( T + 2^ + ^) 3 - 7 2 , (35) 
7 = ^(^ + «-" 2 ) , (36) 
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the system having two excited states. At uj = 0, the lower solution E go , 
corresponding to the ground state \go) = |ll£' 9o 4)o, is 3 



E go = -3G - VG 2 + 85% . (37) 

Using this value, the mean-field parameter 5o = 2G±(go \N+P—\go) can be 
easily calculated, and one obtains 



*o = V 8G 4 " G 2 /8 . (38) 
Differentiating E u with respect to u, the mean value of To can be found 

„(ll£4|r„|ll£4)„ = — ^|j„, co ,«, = J»-w , (39) 

= 4(2G + ^ 
8G<52 - x{2G + x) 2 

The factor J u can be interpreted as the moment of inertia associated with 
the rotations around the Z axis in isospace. It is not defined as positive, 
however for the ground state at uj = (N v = Z v ) it proves to be positive: 

jo = 326 o . ( 41 ) 

If uj ^ 0, but small, the energy E gu can be approximated as: 

tf*,=£flo--y- , (42) 

and the ground state energy £ 9uj , defined by £ 9u = (g w \H\g u ), becomes 
within the above approximations 



£ 9 „=E go + ^_ • (43) 



The parameter uj is calculated from the equations (|27j). (|39j). and for small 
values (iVy — Z v )/2 given to (^ITbl^), it will be 

^W = j^ . (44) 
JO 2J0 ^ ' 



3 For the first excited state Eq. 13311 yields E\ a — —2G, while the corresponding state 
vector CHJ is |1LEi 4} = (|1,-1)|1,1) - |1, 1)|1, -1»/a/2. 
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Replacing uj from (|i3*)l by this value, the energies £ g (A v ; N, Z) of the even- 
even nuclei from the isobaric family having A v = 4 can be expressed as 

£ g (A;N,Z) = E go + {N ~ j p 2 . (45) 

Using these results, some numerical estimates of the constants G4 and 5q for 
the level ld 3 / 2 will be made in the following. In this case, the state |lli? 50 4)o 
having N v = Z v = 2 corresponds to the isosinglet ground state of 36 Ar, 
while |H£s4) w , u = (N v - Z V )/2J$, for N v = 4, Z v = and N v = 0, Z v = 4 
correspond to the T = 2 ground states of 36 S and 36 Ca, respectively. As for 
the members of the same isomultiplet, the ground state energies £ 9 (4;4, 0) 
and £ g (4;0, 4) of 36 S and 36 Ca are the same, higher than the ground-state 
energy £ ff (4;0,0) of 36 Ar. The difference £ g {4; 4,0) - £ g {4; 0,0) is just the 
excitation energy E +. 2 of the state (J 71 "; T) = (0 + ; 2) of 36 Ar at 10.858 MeV. 
So, from a first estimate of J® can be obtained: 

(46) 
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Now, the constant G4 can be calculated from the second-order equation 
given by (|41j) . Following [S] in taking G = 0.6 MeV, two solutions for G4 
are found: G*p = 2.56 MeV, G^ ] = 0.08 MeV. The solution G^ ] is closer 

to the value of the quadrupling interaction constant Gq = 0.09 MeV [S], 

(2) 

suggesting that G4 = G\ is the correct choice, if both terms, Hq and 
come out from the same basic four-particle interaction. For this choice the 
coupling constant of the phenomenological gauge restoring interaction in 
36 Ar has the value 5 = 0.07 MeV. 

In the case of other j = 3/2 valence levels the numerical values of the 
parameters G, So, must be changed because they are dependent on the mass 
number A of the T = and T = 2 isobaric isomultiplets. The average 
will be used for the pairing strength G = 21 /A MeV of the proton and 
neutron constants G p = 23/A MeV, G n = 19/A MeV. The ^-dependence 
of the mean-field parameter 6q is unknown, but it will be obtained from the 
expression of the isorotational energy. The term (A^ — Z) 2 /8J g appearing in 
(|45|) represents the contribution of the nucleons lying on the same degenerate 
j = 3/2 level to the symmetry energy from the Weiszacker mass formula, 
W sym (A, N — Z) = k w (N — Z) 2 , k w = 28.1/ A MeV. This formula describes 
the average symmetry energy of the nuclei, and in particular, for N — Z = 4 
it is expected to approximate the energy E^ + . 2 in the light nuclei having 
N = Z = A/2. A comparison with the experimental data shows that 
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Figure 3. A comparison between the experimental energies E^ + . 2 (* symbols), and 
W sym (A,4) (solid line). 

W sym (A, 4) is systematically greater than E^ + . 2 (figure 3), but except for the 
mass region A < 28 the relative difference W sym (A, 4) / , 2 — 1 decreases 
with A, becoming 1.5 % for 52 Fe. Consequently, for A > 28, W sym (A,4) 
reproduces quite well the value of E^ + . 2 , proving in particular that J° is 

proportional to A. As G ~ A^ 1 , the formula (|41j) also implies an average 
A -1 dependence for 5q, which in the range A > 28 is determined by its value 
at A = 36 to be 5 = 2.7 /A MeV. 

3. The Fermi levels oscillations in superfluid nuclei 

The ground state of the superfluid systems contains correlations between 
the same kind of particles, and is well approximated by the BCS function. 
For a single proton level the pairing Hamiltonian H p and the BCS function 
become: 

H p = eNp - GP+P- (47) 
\BCS) {V)X) = e 2zP +- 2z * p -\0) = e- l ^' 2 \BCS)^ x) (48) 

where A is the Fermi energy and z = pe~ lip . Explicitly, \BCS)^ ^x) 1S 
determined as the ground-state eigenfunction for the Hamiltonian H' = 
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HpL — AiVp, containing the linearised Hamiltonian H p l: 

H pL = eN p - A(P+ + P_) , (49) 

A = G{BCS {0 , X) \P + \BCS {0 , X) ) = — sin4p . (50) 

The term —XN p appearing here can be interpreted as a cranking term restor- 
ing the symmetry broken by H p l. To reach this interpretation, one starts 
from the time-dependent Schrddinger equation idt\^(t)) = H p i\^{t)) writ- 
ten for the Hamiltonian H pL (t) = e - ^*)"^/ 2 H pL e iip ^ N p/ 2 , generated by a 
uniform gauge rotation with the angular velocity = 2 A from H p l- Its so- 
lution |*(t)) = e-^W^/ 2 |* (*)), |*o(*)) = e- iE '^\BCS) {0A) contains the 
previous eigenfunction of H', \BCS)^ 0X ^, whose eigenvalue E' g and parame- 
ter p are determined by 

H'\BCS) {0;X) = E' g \BCS) {0;X) (51) 

to be 



Ef g = -ECl(l-—) , E=^(e-\r + A2 = — , (52) 

cos4p = 2^^ . (53) 

Also, the energy function £ and the mean number of particles M of the 
system are: 

£ = (BCS {0 , X) \H P \BCS {0 , X) ) =en-^-^j- + I ^f + G(Af/Sl), (54) 

M=(BCS (0tX) \N p \BCS {0:X) ) = n-eI p + \I p , (55) 

where I p = 2/G is the "moment of inertia" for the rotation generated by 
N p , and 0(N /Q) = — (e — A — GO/2) 2 /GO is usually neglected, assuming 
O large. A direct calculation shows that similar formulae are obtained for 
the case of 20, levels distributed with a constant density p = 20/ (€2 — ei) 
between the energies e 2 > e\. 

AA c = 0-e^ + A/^ , (57) 
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except for the moment of inertia which becomes I£ = p j coth(2 j pG) |12j . 

To investigate the effects of the interaction —5q{N + P_ + N_P + ) for a 
superfluid system, the microscopic Hamiltonian Hl of 1)20(1 will be treated 
semiclassically 4 , using the variational method for the time-dependent trial 
function |^). If {aj, i = 1, n} are the real time-dependent parameters of \^), 
then the variational equations: 

6 a . J dt{tf\id t - H L \V) = (58) 

give as an approximation of the exact solution of the equation idt \^(t)) = 
H L \$>). Explicitly (JSHJ) becomes: 

H B ii h 3 = q- > H = {*\B L \V) , (59) 
i=i a% 

Bij = 2Im{d j V\d i y) , dj = d/daj , (60) 

describing a classical dynamical system on the trial functions manifold |13| . 
Particularly, choosing j^) as a product of two BCS functions, one for protons 
and another for neutrons, |^) = \BCS p )\BCS n ), this dynamical system 
becomes Hamiltonian. For the one-level case, 

|#) = e 2* P P + -2*;p_ e 2*„JV + -2 2 ;jv_| ^ ^ ZT = p Te -^ ,r = p,n 

and the gauge angles together with the mean number of particles j\f T = 
(^\N T \^) = 20 sin 2 2p T , give the canonical coordinates {Q T , Pr}r=p,n'- 

Q T = (f T , P T = ~y . (61) 

Their time evolution is determined by the Hamilton equations: 



with 



CO 1 

H = 2(e- Y){Pp + Pn) + G(l - ^)(P* + P 2 n ) (63) 



25 JP p P n (O - P p ){0 - P n ) cos(cp p - (p % 



4 This method is expected to yield the average properties of the collective excitation 
produced by the interaction term — 5o(N+P- + N-P+) when fl is large, becoming increas- 
ingly accurate when 8o/G increases (M. Grigorescu, Can. J. Phys. 78 119 (2000)). 
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or equivalently, by the solutions {p T (t), f T (t)}r=p i n of the system 

P P = — sin4/> n sin(v9 p - ip n ) (64) 

Pn = — sin 4p p sm{ipp - ip n ) 

ip p = 2e — G — G(Q — 1) cos Ap p — 5q£1 cot Ap p sin 4/9 n cos(</? p — ip n ) 
(p n = 2e — G — G(Q — 1) cos 4p n — <5o0 cot Ap n sin4/5 p cos(ip p — <p n ) . 

If 5q = 0, then p p = 0, p n = 0, and the number of particles is fixed by the 
initial values p p (0), p n (0), without additional constraints to the variational 
equation (|58j) . When 5q ^ it is convenient to use the " centre of mass" and 
"relative" canonical coordinates (Q,P) and ($,p), respectively, defined by: 

Q = 9e + Qh ; P = P p + P n (65) 



$ = Q p -Q n , p= Pp Pn . (66) 



2 

In these coordinates 7i becomes 

GO 1 P 2 

W(P, p, *) = 2( e - — )P + G{1 - -)(— + 2p 2 ) (67) 



-2^(^-p 2 )[(0-|) 2 -p 2 ]cos<l> 

It is appearing clearly that P = A/2, A = J\f p + N n , is time-independent. 
Moreover, the set M = {(Q, P, <£, p), = 0, p = 0} consists of closed orbits, 
distinguished by A. For each A = 0,4, 8, ...,40 there is an orbit in M 
having = 0, p(t) = 0, P(t) = A/2, and Q(t) = Q + 2At, with 2A = 
2[e — (G + 6 )Sl/2] + [G(l - n- 1 ) + <J ]A/2. If ($, p) are not zero, but small, 
the orbits can be found by expanding TC(P, p, <&) around the point (P, 0, 0). 
Within this approximation Ti becomes: 

H = 2( £ -i^!i)P+[G(l --)+*]— + t p 2 + _ , (68) 

and the equations describing the time evolution of the relative coordinates 
($, p) can be easily integrated to: 

$>(t) = <3?o sin u^i , p{t) = pocos uj v t , (70) 
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p = ^o , u v = VZkC . (71) 

The interpretation of the derivatives Qp/2, Qn/2 as Fermi energies asso- 
ciates the solution given above for the N = Z (p = 0) nuclei, to small 
oscillations in opposite directions of the proton and neutron Fermi levels. 
Pictorially, in the Dasso-Vitturi representation [T3|, they correspond to rel- 
ative angular oscillations of the proton and neutron deformed densities in 
gauge space. Such oscillations, having an isovector character, are allowed 
by the diffuseness of the Fermi surface for the superfluid systems. Their 
quantisation |15j leads to a harmonic oscillator spectrum, whose first ex- 
cited state has the energy ujJ'. For numerical estimates of this energy it is 
convenient to replace the constant k obtained above for a single level by the 
realistic value k = 16k w = 449.6/04 MeV, as given by the symmetry energy. 
The constant C will be calculated for N = Z nuclei with half-filled valence 
shells (sin Ap = 1), taking <5 = 2.7/04 MeV as before, and A = GQ/2. Here 
G = 21 /A MeV, and Q » 0.5(3y4/2) 2 / 3 is half the Fermi level de generacy 
for the harmonic oscillator potential |2j. Using these constants, the energy 
of the isovector monopole vibration becomes uj v = 22.8^4 _1//3 MeV. 

4- Conclusions 

The gauge-restoring interaction —5o(N+P- + iV_P + ) introduced in or- 
der to fix the same proton and neutron Fermi levels in superfluid nuclei was 
proved to be strongly connected with the symmetry energy and with the 
four-particle interaction. For the non-superfluid nuclei, its non-commutation 
with To leads to a kinetic energy term corresponding to collective rotations 
around the Z-axis in isospace (isorotations). In contrast to the purely kine- 
matic symmetry energy determined by the Pauli principle in the Fermi gas 
model ^Hl) the isorotational symmetry energy occurs dynamically from the 
proton-neutron correlations. This result is supported by the experimental 
data on light nuclei, which show large mass differences E +. 2 between the 
even-even members of the T = and T = 2 isomultiplets having the valence 
nucleons on the same degenerate level. An energy Eq+. 2 close to the value 
predicted by the symmetry energy term from the Weiszacker mass formula 
in its simplest form, was taken as an indication of the isorotational origin 
of the T = 2 multiplet. This was the case for the even-even isobars with 
A > 28, and the energy E Q +. 2 corresponding to A = 36 was used to estimate 
the constant 5q as 5o = 2.7 /A MeV. Because the microscopic proton- neutron 

5 When A = 4, Q — 2 this vibration mode correponds to the first excited state |ll_Ei4)o, 
with the excitation energy e x = G + (G 2 + 8<5o) 1//2 - The difference e x — tOv is positive but 
decreases with So, and for So > 1.5G it falls below 10 % of e x - 
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interaction is number conserving, 5o was considered to be a mean-field pa- 
rameter. Such a mean field corresponds to the four-particle interaction ex- 
pressed in (|1U|) . whose strength G4 gives 5o self-consistently. Estimating G4 
for 36 Ar, a value close to the quadrupling constant Gq previously used in 
[S] was obtained. This result shows that the isospin quadrupole interaction 
may represent an alternative to quadrupling for the charge-independent 
treatment of the nuclear correlations. 

In the superfluid nuclei the gauge-restoring interaction leads to isovec- 
tor oscillations of the Fermi levels, superposed over uniform isorotations. 
Their energy oj v = 22.8A~ 1 / 3 MeV was estimated for A > 28 nuclei with 
half-filled valence shells, considering the same constant 5q as for the j = 3/2 
shell. Due to the low value of lo v the Fermi levels oscillations are different 
from the 170 A~ 1/3 MeV resonances presented in 0111 El; but can appear 
as low-lying IVMR in high resolution charge exchange or electron scattering 
experiments. Their occurrence is expected for the superfluid nuclei on the (3 
stability line, the further experimental investigation representing a test for 
the gauge angles dynamics presented above. 
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